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ABSTRACT 

We derive new constraints on the mass of the Milky Way's dark matter 
halo, based on a set of halo stars from SDSS as kinematic tracers. Our sample 
comprises 2401 rigorously selected Blue Horizontal- Branch (BHB) halo stars at 
|z| > 4 kpc, and with distances from the Galactic center up to ~ 60 kpc, with pho- 
tometry and spectra drawn from SDSS DR-6. With distances accurate to ~ 10%, 
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this sample enables construction of the full line-of-sight velocity distribution at 
different Galactocentric radii. To interpret these distributions, we compare them 
to matched mock observations drawn from two different cosmological galaxy for- 
mation simulations designed to resemble the Milky Way, which we presume to 
have an appropriate orbital distribution of halo stars. Specifically, we select sim- 
ulated halo stars in the same volume as the observations, and derive the distribu- 
tions P(Vios/Vcir) of their line-of-sight velocities for different radii, normalized by 
the simulation's local circular velocity. We then determine which value of Vcir(r) 
brings the observed distribution into agreement with the corresponding distribu- 
tions from the simulations. These values are then adopted as observational esti- 
mates for Vcir(r), after a small Jeans Equation correction is made to account for 
slight data/simulation differences in the radial density distribution. This proce- 
dure results in an estimate of the Milky Way's circular velocity curve to ~ 60 kpc, 
which is found to be slightly falling from the adopted value of 220 km s~^ at the 
Sun's location, and implies M(< 60 kpc) = 4.0 ± 0.7 x IO^^Mq. The radial 
dependence of Vcir(r), derived in statistically independent bins, is found to be 
consistent with the expectations from an NFW dark matter halo with the estab- 
lished stellar mass components at its center. If we assume an NFW halo profile 
of characteristic concentration holds, we can use the observations to estimate the 
virial mass of the Milky Way's dark matter halo, Mvir = l-Otoi x IO^^Mq, which 
is lower than many previous estimates. We have checked that the particulars of 
the cosmological simulations are unlikely to introduce systematics larger than the 
statistical uncertainties. This estimate implies that nearly 40% of the baryons 
within the virial radius of the Milky Way's dark matter halo reside in the stellar 
components of our Galaxy. A value for Mvir of only ~ 1 x also (re-)opens 

the question of whether all of the Milky Way's satellite galaxies are on bound 
orbits. 

Subject headings: Cosmology: dark matter — galaxies: individual (Milky Way) 
— Galaxy: halo — stars: horizontal-branch — stars: kinematics 



1. Introduction 

The visible parts of galaxies are, in the current paradigm for galaxy formation, concen- 
trations of baryons at the center of much larger dark matter halos, which have assembled 
through hierarchical merging and gas cooling. Understanding the properties of these dark 
matter host halos, their virial masses, concentrations, and radial mass profiles, vis-a-vis the 
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luminous properties of the main galaxy at their center, is crucial for modeling the dynam- 
ics of the central galaxy, for connecting observations of galaxies to large-scale cosmological 
dark matter simulations, and for understanding what fraction of baryons in the halo ended 
up as stars in the central galaxy. In turn, the extended stellar distributions of galaxies, in 
particular the stellar halos of nearby galaxies, offer some of the best probes to test generic 
predictions about the nature of dark matter mass profiles (Navarro et al. 1996). 

The Milky Way and its surrounding halo are of particular interest, as our internal 
position permits the placement of unique constraints on the Galaxy's stellar mass content, 
on its dark matter profile at large radii, and on the 3-D shape of its dark matter halo. Yet, 
our location within the Galaxy also complicates some measurements, such as the extended 
rotation curve of gas in its disk. As a result, the dark mass profile for the Milky Way between 
~ 10 kpc and ~ 100 kpc and the halo's virial mass have not been previously constrained to 
better than a factor of 2 to 3. In practice, it has proven useful to quantify the halo mass 
profile by either a circular velocity curve, Vcir(r), or by the escape velocity curve, Vcsc(r)- 

In previous work, the most common tools used to estimate the Milky Way halo mass 
are the escape velocity and the velocity dispersion profile of the tracer populations (i.e., 
halo stars, or the Milky Way's satellite galaxies and globular clusters). The escape velocity 
provides constraints on the gravitational potential at the relevant position (Little & Tremaine 
1987; Zaritsky et al. 1989; Kulessa & Lynden-Bell 1992; Kochanek 1996). Recent work 
has shown that the total mass of the halo is around 2 x 10^^ M0. Wilkinson & Evans 
(1999) used the velocities of 27 satellite galaxies and globular clusters to find a halo mass 
of ~ 1.91^7 ^ lO^^ adopting a truncated, fiat rotation curve halo model. Sakamoto, 

Chiba & Beers (2003) used a sample including 11 satellite galaxies, 137 globular clusters, 
and 413 solar neighborhood field horizontal-branch stars, along with a fiat rotation curve 
model, to obtain a total halo mass of 2.5ti'^ ^ 10^^ l-^tar x 10^^ Mq, depending 

on the inclusion (or not) of Leo I. Recently, Smith et al. (2007) estimated a halo mass 
of ~ 1.42^Q54 X 10^^ M0, based on a sample of high- velocity stars from the RAVE survey 
and two published databases, and an adiabatically contracted NEW halo model. Battaglia 
et al. (2005, 2006) used a derived velocity dispersion profile to determine a total mass of 
0.5 ~ 1.5 X 10^^ Mq, with some dependence on the model adopted for the halo profile (see 
also Dehnen, McLaughlin & Sachania 2006). 

In order to improve the precision of the mass estimate of the halo, the fundamental first 
step is to begin with high-quality data (tracers with accurate distances and radial velocities), 
augmented with an efficient method of analysis capable of extracting the maximum amount 
of information. Eor example, Wilkinson & Evans (1999) comment that a data set of ~ 200 
radial velocities of BHB stars reduces the uncertainty in the mass estimate of the halo to 
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~ 20%. This goal has proven elusive, even with data samples twice this size (e.g. Sakamoto 
et al. 2003), due to the use of a relatively nearby sample of tracers. However, we also expect 
a large and distant data set to constrain the gravitational potential at different positions, 
which is ultimately a more reliable probe of the dark matter halo. The planned or already- 
underway kinematically unbiased surveys, such as the European Space Agency's astrometric 
mission Gaia (Turon et al. 2005), the Space Interferometry Mission (SIM; Unwin et al. 
2000), the RAdial Velocity Experiment (RAVE; Steinmetz et al. 2006), the Sloan Digital 
Sky Survey (SDSS; York et al. 2000), and the Sloan Extension for Galactic Understanding 
and Exploration (SEGUE; Newberg et al. 2003, Beers et al. 2004), will finally make it 
feasible to obtain the required large sample of tracers with well-measured parameters. 

Blue Horizontal-Branch (BHB) stars are excellent tracers of Galactic halo dynamics, be- 
cause they are luminous and have a nearly constant absolute magnitude within a restricted 
color range (see, e.g., Sirko et al. 2004 for BHB stars in SDSS). SDSS and SEGUE specifi- 
cally targeted BHB stars for spectroscopy; we employ these stars (and other serendipitously 
discovered BHB stars with available spectra, e.g., mis-identified QSOs) to derive precision 
constraints on the mass of the Milky Way's dark matter halo. 

In §2 we describe the assembly of a particularly conservative (i.e., low contamination) 
selection of BHB stars from the spectra of candidate BHBs in SDSS DR-6. Two different 
cosmological simulations, one of a galaxy resembling the Milky Way with a ~ 2 x 10^^ M0 
halo, and one of a Milky- Way-like galaxy simulation with a ~ 1 x 10^^ Mq halo, are presented 
in §3. This section also describes the estimation of the circular velocity curve as a function 
of Galactocentric radius, Vcir(r), obtained by comparing the observed BHB radial velocities 
to the kinematics of the simulations. In §4 we present the resulting estimates of the Milky 
Way's halo circular velocity curve, and discuss implications for the virial mass of the Milky 
Way's halo. Our conclusions are summarized in §5. 



2. Data 

The SDSS is an imaging and spectroscopic survey covering more than a quarter of the 
sky (York et al. 2000). Although its main focus was (and is) extragalactic science, there 
are still large numbers of stars not only imaged, but also targeted spectroscopically. In 
2005 the project entered a new phase (SDSS-II), in which SEGUE, a sub-survey of SDSS- 
II, addresses fundamental questions about the formation and evolution of our own Galaxy 



(http://www.sdss.org). In addition to further (low-latitude) imaging, SEGUE in particular 
provides for a more systematic and extensive acquisition of spectroscopy for stars in the Milky 
Way, from which stellar parameters and radial velocities can be derived from the application 
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of a well-calibrated and well-tested set of procedures, the SEGUE Stellar Parameter Pipeline 
(SSPP; see Lee et al. 2008a,b; AUende Prieto et al. 2008). 

The Sloan Digital Sky Survey uses a CCD camera (Gunn et al. 1998) on a dedicated 
2.5m telescope (Gunn et al. 2006) at Apache Point Observatory, New Mexico, to obtain 
images in five broad optical bands {ugriz; Fukugita et al. 1996). The survey data-processing 
software measures the properties of each detected object in the imaging data in all five 
bands, and determines and applies both astrometric and photometric calibrations (Lupton 
et al. 2001; Pier et al. 2003; Ivezic et al. 2004). Photometric calibration is provided by 
simultaneous observations with a 20-inch telescope at the same site (Hogg et al. 2001; Smith 
et al. 2002; Stoughton et al. 2002; Tucker et al. 2006). This work is based on stellar spectra 
that are part of the SDSS DR-6 (Adelman-McCarthy et al. 2008). 

The radial velocities are taken from the SSPP, which (primarily) uses matches of the 
observed spectra to a set of 908 ELODIE template spectra, corrected to the heliocentric 
standard-of-rest (HSR) frame. Note that the SSPP automatically corrects for a systematic 
offset in SDSS stellar radial velocities of 7.3 km s~^ (see Lee et al. 2008a,b for derivation 
of this offset). After this correction is applied, the SDSS radial velocities exhibit negligible 
systematic errors, and have an accuracy between 5 and 20 km s~^, depending on the S/N of 
the spectrum and the stellar spectral type. 

2.1. Sample Selection 

We aim to select as "pure" a set of true BHB stars as possible, where the contamination, 
e.g., from halo blue-straggler stars, is minimized, even if this selection procedure results in 
a smaller sample. We have not made any effort to construct a complete sample of BHB 
stars. Therefore, we have not simply followed the DR-6 classification procedures, but have 
employed a very stringent approach combining previously established color cuts with a set of 
Balmer-line profile selection criteria. Our technique is similar to that of Sirko et al. (2004), 
but we use slightly different criteria in the adopted color cuts. Throughout this paper all 
magnitudes are corrected for Galactic extinction, and the colors are corrected for reddening, 
both based on the procedures of Schlegel et al. (1998). 

2.1.1. Color Cuts 

We start our selection of BHB stars by adopting the color cuts for identification of BHB 
candidates used in Yanny et al. (2000): 



- 6 - 



0.8 < u-g < 1.6 
-0.5 < g-r < 0.0 

These color cuts are shown in Figure [U the rectangle is the region of BHB candidates. This 
color cut produces ~ 10, 000 BHB photometric candidate^ with existing spectra, but with a 
considerable contamination by both blue straggler (BS) stars and warm main-sequence stars 
(MS). The subsequent spectroscopic analysis for BHB candidates is aimed at eliminating, or 
at least greatly reducing, contamination from such stars. 

2.1.2. Balmer-Line Profile Cuts 

In the ranges of effective temperature considered herein (roughly 7000 K to 10,000 K), 
BHB stars have lower surface gravities than BS stars and higher temperatures than (old, 
halo population) MS stars. The Balmer-line profiles of warm stars are sensitive to both 
gravity and temperature; their analysis provides a powerful method to select BHB stars 
with confidence. We analyze the line profiles after normalizing the continuum for all stars, 
as illustrated in Figure [21 Then we combine two independent methods to identify non-BHB 
St cirs, cis described below. 

The Dq2 'i^s. fm method (Pier 1983; Sommer-Larsen, Christensen 1986; Arnold & 
Gilmore 1992; Flynn, Sommer-Larsen, & Christensen 1994; Kinman, Suntzeff, & Kraft 1994; 
Wilhelm, Beers, & Gray 1999) discriminates BHB stars from BS stars by determining the 
value of Do. 2, the width of the Balmer line at 20% below the local continuum, and distin- 
guishes BHB stars from MS stars by measuring the value of fm, the flux relative to the 
continuum at the line core (Beers et al. 1992; Sirko et al. 2004). Figure [3] shows how D0.2 
can distinguish a BHB star from a BS star. A plot of D0.2 versus fm of the H^ line for stars 
brighter than g = 18, and passing the initial color cuts, is shown in Figure HI The concen- 
tration of stars centered at (fm, D0.2) = (0.23, 25A) represents the BHB stars; the stars with 
larger D0.2 are BS stars, and the remaining stars are MS stars. Figure H] indicates that the 
sample contamination resulting from the color-selection only is rather severe, on the order 
of 50%. The criteria for confirmation of a BHB star, based on the H^ line analysis, are: 

17 A ^ D0.2 ^ 28.5 A, 0.1 ^ fm ^ 0.3 . 



data file of the full set of candidates is made available in electronic form in the online edition of the 
ApJ, see Table 1 for example. 



- 7- 



The scale width vs. shape method (Clewley et al. 2002) is based on a Sersic profile 
(Sersic 1968) fit to the Balmer hnes: 



where y is the normahzed flux density and Aq is assumed to be the nominal wavelength of the 
Balmer line. To account for small radial velocity corrections and the imperfect normalization 
of spectra, we fit the normalized extracted spectrum to the Sersic profile with five free 
parameters : a, b, c, Aq and n. 



The set of stars that passed the initial color cuts exhibits a bimodal distribution in the 
c^ versus b^ plane (Figure [5]), where 7 refers to the line. This gap allows one to quite 
cleanly separate BHB stars from BS stars, according to: 



As the color coding of the points in Figure \5\ shows, most BHB stars selected by the 
-D0.2 & fm niethod already lie in the appropriate region of the scale width vs. shape method 
applied to H^. The combination of these two stringent criteria indeed appears to eliminate 
most stars that are not bona-fide BHB stars. The combination of Figure H] and Figure [5] 
suggests that the contamination is well below 10%. 

Note that the two spectroscopic criteria are determined based on the spectra of bright 
stars {g ^ 18). Because of the high quality spectra available for such stars, the criteria can 
be identified easily by eye from Figure [Hand Figure [51 For fainter stars {g > 18) we adopt 
the same spectroscopic criteria as for the bright stars. 

There are a total of 2558 starfl that survive the color cuts, and both of the Balmer-line 
profile cuts described above. This sample forms the basis of our remaining analysis. 




(1) 




(2) 



0.75 ^ c^ ^ 1.25 
7.5 A ^ b^ ^ 10.8 - 26.5 (c^ - 1.08)^ . 



data file of the full set of adopted BHB stars is made available in electronic form in the online edition 
of the ApJ, see Table 2 for example. 
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2.2. The Absolute Magnitude of BHB Stars 

Our basic approacli is to identify BHB stars from their spectra (§2.1), then estimate 
their absolute magnitude from photometry alone. BHB stars have similar, but not identi- 
cal, absolute magnitudes, as they are affected slightly by temperature and metallicity (e.g., 
Wilhelm, Beers, & Gray 1999; Sirko et al. 2004). Figure E] shows five theoretical absolute 
magnitudes Mg = 0.60,0.55,0.65,0.70,0.80 in the (u — g, g — r) plane, taken from Sirko et al. 
(2004). For each BHB star in our sample we define the most probable absolute magnitude 
associated with its (u — g, g — r) colors by simply finding the absolute magnitude of thepoint 
on the theoretical track that is closest to the observed star in this color-color spaceo The 
absolute magnitude error of a given star derived from this method is on the order of 0.2 mag 
(Sirko et al. 2004), corresponding to a distance accuracy of 10%; errors in the measured 
photometry are much smaller than this error. The distance from the Sun, and from the 
Galactic center can then be determined from: 

^ = M, + 51ogiorf-5 (3) 

= (R© - dcosb cos if + £ sin^ b + £ cos^ b sin^ I, (4) 

where g is the extinction-corrected magnitude in the g band, is the absolute magnitude 
in the g band, d is the distance to the Sun, r is the distance from the Galactic center, b and 
/ are the Galactic latitude and longitude respectively; we take R©, the distance of the Sun 
from the Galactic center, to be 8.0 kpc. 

2.3. The Spatial, Velocity and Metallicity Distribution of the BHB Star 

Sample 

For ease of the subsequent analysis, we convert the heliocentric radial velocities to the 
Galactic standard of rest (GSR) frame by adopting a value of 220 km s~^ for the Local 
Standard of Rest (Visr) and a Solar motion of (+10.0, + 5.2, + 7.2) km s"^ in (U,V,W), 
which are defined in a right-handed Galactic system with U pointing towards the Galactic 
center, V in the direction of rotation, and W towards the north Galactic pole (Dehnen & 



■^The method described above could be improved upon by fitting a polynomial to the theoretical absolute 
magnitudes and color index (u — g,g — r), and use it instead of picking the "closest" model. However, the 
difference between the two methods is less than the theoretical uncertainties in the derived luminosities. 
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Binney 1998). Hereafter, Vios stands for the radial velocity in the GSR frame (i.e., the radial 
velocity component along the star-Sun direction, corrected for Galactic rotation). If Vhciio 
is the heliocentric radial velocity, then: 

Vios = Vhciio + 10.0 km s""*^ X cos /cos 6 + 
7.2 km s""^ X sin6 + (Visr + 5.2) km s""*^ x sin /cos 6. (5) 

Our sample of 2558 BHB stars may contain some thick-disk stars, with 1 kpc < |z| < 
4 kpc, so we impose an additional geometric constraint |z| > 4 kpc, which reduces the sample 
to 2401 (presumed) halo BHB stars within 60 kpc, and with radial velocity error less than 
30 km s~^. 

The spatial distribution of our 2401 halo BHB stars is shown in Figure [7] and Figure [H 
with all stars located at least 4 kpc from the Galactic disk and at 5 — 60 kpc from the 
Galactic center. Figure M (upper panel) shows the derived [Fe/H] estimates for the BHB 
sample reported by the SSPP (see Lee et al. 2008a,b for details). Most of our BHB stars 
are metal-poor ([Fe/H] ~ —2), as expected for a sample of halo stars. This is further 
testament to the quality of the SDSS spectra and the rigor of the sample selection. The 
observed distribution of line-of-sight velocities is well-fitted by a Gaussian distribution with 
(Tvi^^ = 105 km s~^, as shown in the lower panel of Figure [91 The upper panel of Figure [TO] 
shows the distribution of Vios vs. r for the halo BHB stars. It reveals a nearly equal number 
of stars with positive and negative Vios at a given radius: the BHB population exhibits very 
little net rotation, so we will subsequently analyze only Vios- The lower panel of Figure [TOl 
shows the binned velocity dispersion, cxios, of the sample as a function of radius. It is well 
described by 

'7ios(r) = o"oexp ^ — ^ , with ctq = llll} km s~^ and ro = 3541qq kpc (6) 

where r is the Galactocentric radius. We will use Eqn 6 subsequently for Jeans Equation 
modeling. A maximum likelihood fit to estimate cxo and ro yields (Tq = 111^} km s^^ and ro = 
4271^108 kpC' which is consist with Eqn 6 within 1-a. This indicates that the choice of the 
bins has small effect on the inferred (Tios(r)- We also checked that these values of crios(r) are 
insensitive to small changes in the assumed Visr. 

The radial number density profile of halo stars in the Milky Way in the range ~ 10 — 60 
kpc can be approximated by p ~ ^-3.5 (^g ^1 2007). Accounting for the r^dr volume 

effect, inspection of Figure [TOl and Figure [T2] reveals that the radial distribution of our BHB 
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sample falls off much more rapidly than this, in particular at large distances. This is mostly 
attributable to the SDSS spectroscopic target selection. The chances of a candidate BHB 
star to be targeted for spectroscopy in the course of SDSS and SEGUE is not a simple 
function of its apparent magnitude, but depends on many (often operational) factors, which 
have also evolved over the course of the survey. The net result is that only under very 
favorable circumstances are distant (hence faint) BHBs targeted, and have spectra obtained 
that are of sufficient S/N to pass our quality criteria. As a result, the radial distribution of 
our sample falls off much more steeply than the parent population of halo BHB stars. We 
account for these effects in the subsequent analysis, as described below. 

Overall, our present sample of distant halo stars with available kinematics is nearly an 
order of magnitude larger than that of Battaglia et al. (2005), which has 240 stars. Because 
our tracers are BHB stars, their distances are also known more accurately. However, the 
Battaglia et al. sample does extend to larger radii, up to ~ 100 kpc. 

3. Modeling the BHB Kinematics 

We now describe our approach to convert the Vios(r) measurements (as a function of 
Galactocentric radius, r) for the ~ 2400 BHB stars into estimates of Vcir(r) of the Milky Way 
halo, and ultimately to estimate M^ii. for the Milky Way's halo. With the full line-of-sight 
velocity data set, sampling radii of 5 — 60 kpc, we can obtain both the velocity dispersion, 
(Tios, and identify a set of exceptionally high- velocity stars at various radii that might be 
suitable to estimate Vesc(r)- Some previous work (e.g., Battaglia et al 2005) considered only 
velocity dispersions, while others focused only on high-velocity stars in the Solar neighbor- 
hood (Sakamoto et al. 2003; Smith et al. 2007). We consider the full velocity distribution; 
as we shall see later, it is close to a Gaussian and hence most of the pertinent information 
is contained in aios- 

To link the observables, Vios and r, to Vcir(r), we must not only account for the particular 
survey volume, but also need to make at least an implicit assumption about the nature of the 
halo-star distribution function, in particular its (an-)isotropy. In contrast to, e.g., Battaglia 
et al. (2005), we restrict ourselves not only to Jeans Equation modeling, but choose to 
account for these issues also by investigation of a comparison with cosmologically motivated 
galaxy simulations, to make "mock observations" within these simulations, and then match 
test results to the observations. One cannot expect the halo stars in the simulations to 
have exactly the sample density profile as the actual stars in the halo of the Milky Way; we 
account for such differences as described below. 
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3.1. Halo Star Kinematics in Simulated Galaxies 

Based on prior estimates of the halo mass, we chose two SPH simulations corresponding 
to the formation of Milky Way-like haloes, from which "pseudo-observations" are constructed 
for comparison with the data. Both halos were picked from low-resolution cosmological dark 
matter simulations and were re-simulated at higher resolution, including the effects of gas 
dynamics, and star formation and stellar feedback. 

Simulation 1 was run using the smoothed particle hydrodynamics code GRAPESPH 
(Steinmetz 1996), and included a moderately efficient stellar feedback model in which the 
supernova energy was added to the thermal energy of the gas, as described in Abadi et al. 
(2003). This simulated halo has a present day virial mass of about ~ 8.6 x 10^^ Mq and a 
virial radius of 206 kpc. 

Simulation 11 was run using the smoothed particle hydrodynamics code GADGET-2 
(Springel 2005), assuming the standard star-formation prescription of Springel et al. (2003), 
but without stellar feedback. This simulated halo has a present-day virial radius of rvir = 345 
kpc, and a virial mass of Mvir = 2.1 x IO^^Mq (these values differ from Naab et al. 2007, as 
we use a different virial contrast in this paper). 

For further details of the simulations we refer the interested reader to Abadi et al. (2003, 
Simulation I) and to Naab et al. (2007, Simulation 11), and references therein. 

For each particle in these simulations we have the 3-D positions and 3-D velocities, the 
circular velocity, Vdr = a/ id^/di, and the escape velocity, Vcsc, to the virial radius of the 
simulation at each given position. Because SDSS/SEGUE only observed ~ vr steradians, not 
the entire sky (see Figure [TTI) . in order to compare the observations to the simulations we 
must select simulated star particles in the same region as the SDSS footprint. 

As a first step, we specified "Galactic coordinates" in the simulation by defining the 
"Galactic Plane" by the net angular momentum of all (stellar) particles within 10 kpc of 
the center. In this coordinate system one has (by definition) Ltot = \^tot\Gz, and the Sun is 
at (x0,y0,O.O) kpc (R© = y/x^ -|- y| = 8.0 kpc). We can then calculate the Galactocentric 
radial velocity as seen from the "Sun" for each simulation particle as: 



= v/(x - Xq)2 + (y - yQ)2 + z2 
Vx(x-Xq) +Vy(y-y0) +v,z 

Vios - 2 ^ ' ' 

Since the simulated "observer" is at rest in the Galactocentric coordinate system, no as- 
sumptions about the LSR need be made in the simulations. 
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To create "pseudo-observations", we first remove the satellites in the simulation and 
then assume 10 positions for the Sun that all have R = 8 kpc and z = kpc, but different 
azimuthal angles, 0. The corresponding / and b in this coordinate system can then be used 
to select the star particles in the same region of the simulation as our observed BHB sample. 
The effective flux limit of the BHB sample, g < 20, implies a maximum distance from the 
Sun of about 76 kpc, a distance which we impose as a selection limit in the simulations. To 
select "halo" stars, we also impose that the stars have to be at least 4 kpc above or below 
the disk plane. Finally, we average these 10 samples of simulated halo stars. Figure [TT] shows 
the distribution of these selected simulated stars and the observed halo BHB stars in I and 
b. 

This procedure results in a sample of simulated halo stars (each position produces an 
"observational data set") with Galactocentric radial velocities, Galactocentric radii, escape 
velocities, and circular velocities, whose distribution is shown in Figure [T2l This figure makes 
it clear that, even in a large sample, the Galactocentric radial velocity rarely approaches the 
escape velocity (one obvious reason being that we measure only the projected component of 
the space velocity). 

3.2. Estimating Vcir(r) from the data 

We analyze the implications of the observed BHB kinematics for the Milky Way halo's 
mass distribution in two steps. We first estimate Vcir(r) from the data/ simulation comparison 
in a set of statistically independent radial bins, effectively constructing a circular velocity 
curve extending to 60 kpc. We then fit the circular velocity curve with NFW halo (and 
bulge+disk) models, resulting in estimates of Mvir. It should be noted that the Milky Way 
halo's presumed virial radius extends about a factor of four beyond the most distant BHB 
stars in our observed sample of tracers. 

For the data/model comparison we construct the distributions, P(Vios/Vcir(r)), for the 
simulations and the data respectively. For the simulations, we use the procedure described 
in the above section to obtain the Vios and Vcir at each particle position. We compare those 
Psim(Vios/Vcir(r)) distributions to analogous ones constructed from the data for a sequence 
of trial values Vcir(r). As the best observational estimate of Vcir(r), we then take that value 
for which the probability that Fobs (Vios /Vcir) and Psim(Vios/Vcir) were drawn from the same 
distribution is maximal. As we have no a-priori functional form for these distributions, we 
define this best match as that which maximizes the probability in a two-sided K-S test. To 
define confidence limits on the Vcir(r) estimate, we repeat this procedure with bootstrapped 
versions of Fobs ( Vios/ Vdr), and take as 6^^.^ the variance of the resulting V^r distribution. 
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We divide the Galactocentric radius (r) into 10 bins and apply this method, which results 
in 10 (statistically) independent radial estimates of Vcir(r). Specifically, the bins we adopted 
were 5-10 kpc, 10 - 15 kpc, 15 - 20 kpc, 20 - 25 kpc, 25 - 30 kpc, 30 - 35 kpc, 35 - 40 kpc, 
40 — 45 kpc, 45 — 50 kpc, and 50 — 60 kpc (all distances are from the Galactic Center, not 
from the Sun). The best matched distributions of P(Vios/Vcir) are illustrated in Figure [T3] 
and Figure [m In these figures, the red line shows P(Vios/Vcir) of the simulations, while the 
black line shows that of the data. Overall, the velocity distributions agree well with one 
another. 

To estimate the statistical uncertainties on these values, we use bootstrap resampling 
on our BHB sample (typically 100 times), repeat the above procedure, and indicate the 
resulting 68% (1-cr) confidence region as error bars. 

In a given gravitational potential, more centrally concentrated kinematic tracer popula- 
tions will exhibit smaller velocity dispersions than more extended ones. In the simple case of 
a spherical potential, with tracers of p ~ r"''' and an isotropic velocity dispersion that varies 
only slowly with radius, the Jeans Equation, for a given Vcir(r), yields the relative velocity 
dispersions of the two populations, 01/02 = ^72/71, where the indices refer to two different 
hypothetical tracer populations. 

In the radial range of 10-60 kpc, the density profile of halo stars in Simulation II is 
approximately p ~ r~^'^, while that of Simulation I is p ~ ^--3.7 This, however, should 
not be compared to the radial distribution of the stars for which we actually have velocities, 
but to the radial profile of halo stars from which they were drawn (see §2.3). This requires 
the reasonable assumption that the measured velocities are uncorrelated with the spectral 
targeting - even a more complete sampling of BHB stars at a given radius would have yielded 
the same velocity distribution. As an estimate for the actual density profile for the Milky 
Way's stellar halo at 10-60 kpc, we take the estimates (for main sequence turn-off stars) 
of p ~ r~'^-^ (Bell et al. 2007), based on SDSS. Using the above correction based on the 
Jeans Equation, we must subsequently revise the derived velocity scales for the Milky Way 
halo upward by ^^3.5/2.9 = 1.1 for Simulation II and downward by ^^3.5/3.7 = 0.97 for 
Simulation I. 

Matching P(Vios/Vcir(r)), and applying the above correction, we obtain the estimates 
summarized in Figure [T5l and in Table 3. The filled circles in this Figure refiect the Vcir(r) 
estimates based on Simulation I, while the filled squares stand for the Vcir(r) estimates based 
on Simulation II. As mentioned before, the error bars are from bootstrapping. 

We have explored whether the derived estimates of Vcir(r) depend on our adopted value 
for the distance of the Sun from the Galactic center, R©, or on the adopted local rotation 
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velocity, Visr- If we vary Rq adopted in the observation and simulations (see Eqn 4 and Eqn 
7), from 7.5 kpc to 8.5 kpc, Vcir(r) changes by only about 1%. Our results for the derived 
Vcir(r) are also insensitive to the choice of the Visr adopted in the observation (see Eqn 5). 
Taking Visr as 200 km s~^ or 240 km s~^ changes the estimated Vcir(r) by only 3%. 

For reference, we show how these estimates of Vcir(r) compare to those derived from the 
Jeans Equation and the fit to o"ios(r) shown in Figure [TOl From the Jeans Equation, Vcir(r) 
can be estimated from the velocity dispersion, cTj. (Binney & Tremaine 1987) as follows. 



r dja^p) 
p dr 



2Pa^ = Vlir) 



with 



(9) 



where o"r(r) and o"t(r) are the radial and tangential velocity dispersions, respectively, in 
spherical coordinates, and p(r) is the stellar density. 

The distribution of the halo stars in the simulations are anisotropic, with (3 = 0.37, 
and the simulations exhibit o"ios(r) ~ cri.(r) for this particular survey volume. Taking from 
Figure [To] the fit to crios(r) of the BHB stars (Eqn 6), and assuming (Tios(r) ~ o"r(i') for BHB 
stars, we can derive two circular velocity curves, for (3 = 0.37 (anisotropic), and (3 = 
(isotropic), respectively, by adopting p(r) ~ r~^-^ (see Figure [15]). They both agree well 
with the simulation-based estimates, and are not used any more in the following analysis. 



4. Results 



Figure [TBI and Figure [T7I present the main result of our analysis, Vcir(r), an estimate of 
the circular velocity curve from ~ 10 — 60 kpc. This represents the first time that the circular 
curve for the Galaxy has been estimated to such large distances at this accuracy. Note that 
at small radii this estimate, though derived from halo stars, agrees well (within ~ 10%) with 
established determinations at the Solar radius (~ 220 km s~^). Beyond the Solar radius, the 
circular velocity curve appears to be gently falling to 175 km s~^ at ~ 60 kpc. Note also that 
the circular velocity curve is a conceptually more robust estimate than Vesc, which depends 
more sensitively than Vcir on $(r) at radii beyond the measurements. 

Using the functional form for Vcir(r) expected for an NEW halo and the stellar compo- 
nent (see below) as a means to interpolate the individual circular velocity curve estimates, 
one obtains Vch-(60 kpc) = 170 ± 15 km s-\ or M(< 60 kpc) = 4.0 ± 0.7 x 10"Mq. This is 
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the largest radius for which the data directly constrain Vcir(r) or M(< r). Yet this radius is 
only one-fourth of the expected virial radius of the Milky Way's halo. We therefore proceed 
with a separate step, to use these Vcir(r) estimates to constrain parameterized models for 
the overall dark matter halo. We assume that the Galactic potential is represented by three 
components, a spherical Hernquist (1990) bulge, an exponential disk for the stellar compo- 
nents, and describe the halo by an NFW profile (Navarro et al. 1996). The total potential 
can then be simply expressed as 

<l>tot(r) = $disk(r) + $buigc(r) + $NFw(r), (10) 
with an assumed potential, presumed to be spherically symmetric, for the disk and bulge of 

GMdisk(l-e-^) 

$disk(r) = , (11) 

r 

*bulge(r) = -^Mi^, (12) 
r + Co 

where Mbuige = 1-5 x 10^° Mq, cq = 0.6 kpc, Mdisk = 5 x 10^° M©, and b = 4 kpc (similar to 
Smith et al. 2007). The radial potential for a spherical NFW density profile can be expressed 

as 

W)^-i^ln(l + ^). (13) 

where c is a concentration parameter, defined as the ratio of the virial radius to the scale 
radius. For standard ACDM cosmogonies we do not attempt to constrain halo flattening. 
The parameter ps is a characteristic density given by 

~ 3 ln(l + c)-c/(l + c)' ^ ' 

where pcr = SH'^/SirG is the critical density of the universe, is the contribution of matter 
to the critical density, and 5th is the critical overdensity at virialization. The virial mass can 
then be determined from the virial radius, using 

Mvir = — Pcr^m^thfvir- (15) 

For our analysis we adopt flja = 0.3, 5th = 340, and Ho = 65 km s~^ Mpc~^. Given recent 
discussions (and doubts raised) of whether the baryons modify the dark matter profile, as 
expected from "adiabatic contraction" (Button et al. 2007), we consider both an unaltered 
and an adiabatically contracted NFW profile in the fit of $tot- 
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By fitting the observed Vcir(r) with A/rd$/dr from $tot(r) shown as Eqn 10 we can 
constrain the halo mass of the Milky Way. In this fit, we simply adopt an unaltered NFW 
profile and a present-day relation between the mean value of c and Mvir, 



logio c = 1.075 - 0.12(logio M,ir - 12). (16) 

This relation is accurate over the range 11 < logMyij. < 13, and is based on the model of 
Maccio et al. (2007) with = 0.3, Q\ = 0.7, as = 0.9, and Uj = 1.0. Therefore, the Mvir is 
derived as a 1-parameter fit (fit only Mvir, presuming c(Mvir))- The results are summarized 
in Figure [T6l 

Specifically, for the circular velocity estimates resulting from Simulation I, we find 

Mvir = 0.91^°;?^ X 10^2 Mo, 

with rvir = 267^19 kpc and c = 12.0^0:3. 

For the circular velocity estimates based on Simulation II we find 

Mvir = 0.82t^oil X 10'' Mq, 
with rvir = 2581^; kpc and c = 12.2l°j. 

Note that the error bars of Mvir, fvir, and c are determined from 1-a confidence intervals 
in a Chi-Square test. For the two cases above we have adopted an unaltered NFW profile 
and an average relation between c and Mvir- If we fit an adiabatically contracted NFW 
profile (using the prescription of Blumenthal et al. 1986 and Mo et al. 1998) and the same 
disk and bulge as in Eqn 11 and Eqn 12, taking the concentration parameter (c) and virial 
mass (Mvir) of the NFW profile as independent parameters (i.e., we do not require that they 
follow the relationship in Eqn 16), the Mvir can be derived as a 2-parameter fit (Mvir and c), 
as shown in Figure [T71 

For the circular velocity estimates resulting from Simulation I 

Mvir = l.Otoi X IO^^Mq, 

with rvir = 275120 kpc and c = 6.6;||. 

The virial mass calculated by the circular velocity estimates based on Simulation II is 

Mvir = 1.2ll|]:^[J X lO^'Mg, 

with rvir = 293^26 kpc and c = A.Stli- 
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The error bars on M^ir and rvir are determined for 1-a confidence intervals in a Cfii- 
Square test, fixing the best-fit value of c, while the error bars on c are determined for 1-a 
confidence intervals in a Chi-Square test, fixing the best-fit value of Mvir- 

Note that although contracted and uncontracted halo fits differ quite strongly in their 
(initial) concentration, reassuringly the Mvir estimates remain relatively unaffected. The 
lower concentrations for the contracted halo fits are still reasonally consistent with the con- 
centration scatter expected from cosmological simulations (e.g., NFW 96). 

The results here show that the Milky Way's circular velocity curve must be gently falling 
to distances of 60 kpc from its value of ~ 220 km/s at the Solar radius; the null hypothesis 
Vcir(r) has a constant value of 220 km is rejected by our fits with very high statistical 
significance (at a level of 0.01). 

A direct comparison with earlier work, at the data or crios(r) level, is not straightforward 
to carry out, as each sample has distinct selection effects, such as differing radial distributions. 
Our estimate of Mvir, taken to be the average of the four estimates derived in this Section, 
1.01q'2 ^ IO^^Mq, falls well within the (considerably larger) range of values estimated by 
Battaglia et al. (2005, 2006): 0.5 ~ 1.5 x 10^^ M©. This value is also reasonably consistent 
with the recent estimate by Li & White (2008), based on Local Group dynamics. In general, 
however, the new mass estimate presented here lies at the lower limit of most previous 
estimates. 

The method to estimate Vcir(r) can also be used to derive the escape velocity curve 
Vesc(r). We have carried out an analogous procedure to that described in §3.2, using 
P(Vi.o.s/Vesc), and found lower hmits on Mvir > 0.5 x which do not appear to 

be as stringent as the Vcir(r) estimates. 

5. Summary and Conclusion 

We have constrained the mass distribution of the Milky Way's dark matter halo, by 
analyzing the kinematics of nearly 2400 BHB stars drawn from SDSS DR-6, which reach 
to Galactocentric distances of ~ 60 kpc. To obtain a "clean sample" of BHB stars, we 
have re-analyzed all candidate BHB spectra, following the prescription by Sirko et al. 2004, 
which should result in a contamination fraction (mostly by BS stars) of well below 10%. The 
metallicity distribution, centered on [Fe/H] ~ —2 dex, confirms that most sample members 
must be halo stars. The distances to the BHB stars are known to ~ 10%. 

To account for the complex survey geometry and for plausible orbital distributions of 
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our sample of BHB stars, we have compared the observed radial velocities to analogous 
quantities drawn from the "star particles" in galaxy formation simulations that resemble the 
Milky Way. In particular, we have placed a virtual "observer" at 8.0 kpc from the simulation 
centers, looking in the actual SDSS directions and sampling radial velocities for stars that 
are more than 4 kpc above and below the disk plane. We then explored to what mass scale, 
or Vcir(r), the simulations need to be scaled to in order to match the observed line-of-sight 
velocity distribution in a set of radial bins. In this analysis, we adjusted this scaling using 
the Jeans Equation, to account for the slight difference in the radial profile of observed halo 
stars (p ~ simulated halo stars (p ~ ^^^^ ^^lis radial range. 

This procedure results in direct estimates of Vcir(r) from ~ 10 kpc to ~ 60 kpc, the best 
such estimates to date over this range. The circular velocity estimate varies slightly with 
radius, dropping from ~ 220 km s~^ at 10 kpc to ~ 170 km s~^ in the most distant two bins. 
Applying this procedure to two independent cosmological simulations (Simulation I and Sim- 
ulation II, respectively) results in consistent estimates of Vcir(r). The mass enclosed within 
60 kpc, constrained quite directly by the data, is found to be 4.0 ±0.7 x lO^^M©. As a result 
of the much more extensive data set provided by SDSS, the uncertainties on this estimate 
are substantially lower than those obtained by previous comparable work. A simple, Jeans- 
Equation based modeling approach, assuming (an-)isotropies of either (3 = 0.37 or /? = 
(as found for the halo stars in the cosmological simulations) yields results that are consistent 
with these values. 

Although each of the Vcir(r) points were estimated independently, the implied overall 
profile is consistent with both the mass profile in the simulations and with a parameterized 
mass model that combines a fixed disk and bulge model with an NEW dark matter halo, 
whose concentration c corresponds to the expected mean value for its virial mass. We have 
also explored halo fits with the concentration c as a free parameter, as well as halo profiles 
that have been modified by adiabatic contraction. We find that our data cannot discriminate 
well whether adiabatic contraction occurs or not - an uncontracted halo of higher than 
average concentration and a contracted halo of (initially) low concentration fit comparably 
well. The resulting virial masses, are consistent for both fitting 

approaches. We have also checked that these results are quite robust with respect to distance 
errors, modest sample contamination (< 10%), and the choice of a different, independent 
galaxy formation simulation. 

The estimate of Myir, which does imply an extrapolation from r^ax = 60 kpc to the 
virial radius of ~ 250 kpc, is consistent with a recent estimate from a much smaller sample 
of halo stars (Battaglia et al. 2006), but it is lower than previous estimates that also rely on 
the kinematics of satellite galaxies (e.g., Kochanek et al. 1996). However, recent results on 
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the LMC (Kallivayalil et al. 2006; Besla et al. 2007) indicate that not even the Magellanic 
Clouds may have been bound to the Milky Way for long, posing a potential conceptual 
problems for the use of satellite dynamics. It should be interesting to explore how our new 
constraint modifies the Local Group timing argument and its implication for M31's halo 
mass (Li & White 2008). 

The estimate of M^ir ~ IO^^Mq, together with an estimated total cold baryonic mass of 
6.5 X 10^°Mq and a global baryon mass fraction of 0.17, implies that nearly 40% of all baryons 
within the virial radius have cooled to form the Milky Way's stars and (cold) gas, consistent 
with recent estimates for galaxies of that mass scale, based on statistical arguments (van 
den Bosh et al. 2007; Gnedin et al. 2007). 

We note that parts of our analysis have been performed under the assumption that the 
stellar halo of the Galaxy is considered as a single relaxed population, or one that matches the 
simulations. Our data on the overall dynamics are consistent with that hypothesis. Recent 
evidence from Carollo et al. (2007) and Miceli et al. (2007) suggest that the halo may 
well be more complex, and comprise two distinct populations of inner- and outer-halo stars, 
with (slightly) different net rotational velocities and spatial profiles. We defer all analysis of 
kinematic and spatial sub-structure in this particular sample to a future paper. 

The SDSS and SEGUE surveys have shown in this context that they can provide un- 
paralleled sets of kinematic tracers for the Milky Way, enabling a direct "circular velocity 
curve" estimate of the Milky Way extending to 60 kpc. Once the full set of spectroscopy 
from SEGUE is available, a much larger set of stars for such an analysis should be available. 
The proposed extension of SDSS, known as SDSS-III (and using a more sensitive, 1000 fiber 
spectrograph), will provide the opportunity to obtain higher quality spectra for fainter, more 
distant BHB stars, thus extending the reach of our analysis to over 100 kpc. 
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Fig. 1. — SDSS color-color diagram showing all spectroscopically targeted objects that were 
subsequently confirmed as stars. The large Balmer jump of A-type stars places them in the 
region, where our "color cut" selection box is drawn. This color selection approach follows 
Yanny et al. (2000). 
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Fig. 2. — Spectrum of a typical (high-S/N) BHB star (upper panel), and the ~ region 
of the same star with the continuum divided out (lower panel). The (fm,Do.2) parameters 
that are used for the sample selection are labelled for both hues (see §2.1.2 for discussion). 
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Fig. 3. — Normalized spectrum of a BHB star (solid line) and a BS star (dotted line), of 
similar effective temperatures, in the ~ H5 region. Although subtle, one can notice that 
the BS stars's Balmer lines are wider at 20% below the local continuum than the BHB star. 
These effects arise due to the higher gravity of the BS star. 
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Fig. 4. — The parameters fm and D0.2, as determined from the hne, for stars brighter than 
g = 18 and passing the color cuts shown in Figure [H The trail of stars with fm ^ 0.30 are 
too cool to be BHB stars, while the concentration of stars with D0.2 ^ 28.5 A is due to blue 
stragglers with higher surface gravity. The stars that lie well outside the main locus are the 
result of poor parameter determinations due to missing spectroscopic data at the location 
of the line. The region enclosed by the box is used as the BHB selection criterion for the 
H^, D0.2 & fm method. 
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Fig. 5. — Green dots are BS stars and blue dots are BHB stars identified by the D0.2 & fm 
methods. They can be distinguished clearly by the "gap". The closed region indicates the 
scale width-shape criteria that selects BHB stars. 
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Fig. 6. — Color-color (u — g vs. g — r) diagram of our BHB stars. The five black dots, starting 
from the right to the left, represent the model colors for BHB stars of absolute magnitudes 
Mg = 0.60, 0.55, 0.65, 0.70, 0.80. 
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Fig. 7. — The spatial distribution of the BHB stars in our sample, projected in the z-x (upper 
panel) and z-y plane (lower panel), respectively. The large filled circle is the Sun; the two 
lines are the planes 4 kpc above and below the Galactic disk plane. 
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Fig. 8. — The spatial distribution of the sample BHB stars in the x-y plane (upper panel) 
and z-r plane (lower panel). The coordinate origin is the location of the Galactic center. 
The two lines in the z-r plot are the planes 4 kpc above and below the Galactic disk plane. 
The short dashed vertical lines in the lower panel show the bin boundaries chosen for the 
subsequent analysis. 
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Fig. 9. — (Upper panel) The distribution of metallicities, [Fe/H], as a function of apparent 
magnitude, for the entire sample of halo BHB stars. (Lower panel) The distribution of line- 
of-sight velocities, corrected to the GSR, for the entire sample of BHB stars. A Gaussian of 
width a = 105 km s~^ centered on the local standard of rest is shown for reference. 
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Fig. 10. — (Upper panel) The distribution of Vios as a function of Galactocentric distance, 
r, for the entire sample of halo BHB stars. (Lower panel) The velocity dispersion, cxios, as a 
function of Galactocentric distance. A best fit exponentially falling relationship is plotted. 
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Fig. 11. — The Galactic sky coverage of the observed BHB stars (red dots) and selected 
simulated stars (black dots), drawn from Simulation I. 



- 35 - 




Fig. 12. — The Galactocentric radial velocity, escape velocity and circular velocity distribu- 
tions of the stars in Simulation I (upper panel) and in Simulation II (lower panel) (see §3); 
the simulations are "viewed" from the position of the Sun to lie within the SDSS DR-6 foot- 
print. The solid hne delineates the predicted escape velocity, while the dashed line indicates 
the predicted circular velocity. The dots represent the radial velocities. 
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Fig. 13. — Comparison of the Galactocentric line-of-sight velocity distribution, P(Vios/Vcir) 
between the halo star particles in the Simulation I and the observations, shown here for all 
bins. The dashed line is P(Vios,sim/Vcir), and the solid line is P(Vios,obs/Vcir), after finding 
the best matching velocity scaling, Vch, listed in Table 3; the best match in this context is 
obtained from a K-S test. 
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Fig. 14. — Comparison of the Galactocentric line-of-sight velocity distribution, P(Vios/Vcir) 
between the halo star particles in the Simulation II and the observations, shown here for all 
bins. The dashed line is P(Vios,sim/Vcir), and the solid line is P(Vios,obs/Vcir), after finding 
the best matching velocity scaling, Vch, listed in Table 3; the best match in this context is 
obtained from a K-S test. 
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Fig. 15. — The distribution of circular velocity estimates, Vdr, for two different simulated 
galaxies. The filled circles are the Vcir estimates for the observed halo BHB stars based on 
Simulation I and the filled squares are the Vcir estimates based on Simulation II. The two 
lines are the circular velocity curve estimates derived from the velocity dispersion profile 
(Figure fTOl) and the Jeans Equation with (3 = 0.37 and (3 = 0. 
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Fig. 16. — Circular curve estimates matched by a combination of a stellar bulge, disk and an 
an unaltered NFW dark matter profile. The solid line is the best-fit circular velocity curve 
to the Vcir(r) estimates, while the large symbols in the two plots are the Vcir(r) estimates. 
Contributions of the adopted model components (i.e. disk, bulge, and halo) and the circular 
velocity curves based on the Jeans Equation are plotted in different linestyles. Estimates of 
virial mass, Mvir, virial radius, rvir and concentration parameter, c, are labeled on the plots. 
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Fig. 17. — As in Figure [161 but here the circular velocity curves were derived under the 
assumption of a contracted NFW profile. The solid line is the best-fit circular velocity curve 
to the Vcir(r) estimates, while the large symbols in the two plots are the Vcir(r) estimates. 
Contributions of the adopted model components (i.e. disk, bulge, and halo) and the circular 
velocity curves based on the Jeans Equation are plotted in different linestyles. Estimates of 
virial mass, Mvir, virial radius, rvir and concentration parameter, c are labeled on the plots. 
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Table 1. List of 10224 BHB candidates selected by color cut 



SpName 


Name 


RA 


Dec 


1 


b 


g 


g - r 


u 


- g 


uerr 


gerr 


re 








degree 


degree 


degree 


degree 


mag 


mag 


mag 


mag 


mag 


m 


51602-0266- 


■225 


SDSSJ094218.23-002519.7 


145.575943 


-0.422125 


236.196579 


36.896130 


15.46 


-0.10 


1 


.20 


0.02 


0.02 


0. 


51602-0266- 


■397 


SDSSJ094138.17-f 001821. 5 


145.409058 


0.305967 


235.328827 


37.175690 


18.48 


-0.11 


1, 


.22 


0.04 


0.02 


0. 


51602-0266- 


■634 


SDSSJ094840.23-I-002818.0 


147.167633 


0.471673 


236.434402 


38.711887 


18.01 


-0.02 


1 


.02 


0.03 


0.02 


0. 


51609-0292- 


■102 


SDSSJ125223.54-003708.2 


193.098083 


-0.618937 


303.444397 


62.251862 


14.29 


-0.24 


1, 


.27 


0.02 


0.02 


0. 


51609-0292- 


■155 


SDSSJ125050.87-000806.1 


192.711960 


-0.135032 


302.609894 


62.736347 


17.84 


-0.00 


1, 


.06 


0.02 


0.02 


0. 


51609-0292- 


■232 


SDSSJ124759.81-000456.2 


191.999207 


-0.082266 


301.051117 


62.776917 


18.16 


-0.04 


1 


.05 


0.03 


0.02 


0. 


51609-0292- 


■269 


SDSSJ124721. 11-002931. 5 


191.837952 


-0.492089 


300.729431 


62.362186 


19.17 


-0.09 


1, 


.12 


0.05 


0.03 


0. 


51609-0292- 


■329 


SDSSJ124641. 66-1-003751. 2 


191.673569 


0.630884 


300.275543 


63.478165 


17.15 


-0.10 


1 


.23 


0.02 


0.01 


0. 


51609-0292- 


■351 


SDSSJ124449.35-I-002157.4 


191.205612 


0.365958 


299.263092 


63.190571 


17.64 


-0.26 


1 


.13 


0.02 


0.01 


0. 


51609-0292- 


■367 


SDSSJ124805.12-|~010113.5 


192.021347 


1.020428 


301.028198 


63.879784 


16.38 


-0.17 


1 


.22 


0.03 


0.02 


0. 



Note. — The first two columns arc object names and the next four columns contains the astrometry (ra, dec, 1, b) for each object. The magnit 
velocities and errors are listed next. The next four columns are the lincwidth parameters from the Balmcr lines. The last column is the classificatic 
edition of the Journal. The printed edition contains only a sample. 



Table 2. List of 2558 BHB stars selected from SDSS DR6 



SpName 


Name 




RA 


Dec 




1 




b 


g 


u 


~g 


g - 


- r 




Do. 2, (5 










degree 


degree 


degree 


degree 


mag 


mag 


mag 


mag 


A 


51602- 


■0266- 


■125 


SDSSJ094317.57-011021.2 


145, 


.823196 


-1.172550 


237, 


.138931 


36. 


,662491 


16.48 


1, 


,17 


-0, 


,22 


0.55 


26.60 


51602- 


■0266- 


■397 


SDSSJ094138.17-f 001821. 5 


145, 


.409058 


0.305967 


235, 


.328827 


37, 


,175690 


18.48 


1, 


.22 


-0, 


,11 


0.55 


27.78 


51602- 


■0266- 


■634 


SDSSJ094840.23-I-002818.0 


147, 


.167633 


0.471673 


236, 


.434402 


38, 


,711887 


18.01 


1, 


.02 


-0, 


,02 


0.60 


25.80 


51609- 


■0292- 


■102 


SDSSJ125223.54-003708.2 


193, 


.098083 


-0.618937 


303, 


.444397 


62. 


,251862 


14.29 


1, 


,27 


-0, 


,24 


0.55 


27.00 


51609- 


■0292- 


■329 


SDSSJ124641. 66-1-003751. 2 


191, 


.673569 


0.630884 


300, 


.275543 


63, 


,478165 


17.15 


1, 


.23 


-0, 


,10 


0.55 


25.03 


51609- 


■0292- 


■351 


SDSSJ124449.35-I~002157.4 


191, 


.205612 


0.365958 


299, 


.263092 


63. 


190571 


17.64 


1, 


,13 


-0, 


,26 


0.60 


24.54 


51609- 


■0292- 


■582 


SDSSJ125254.01-I-002903.2 


193, 


.225037 


0.484220 


303, 


.747009 


63. 


,353657 


16.22 


1, 


,23 


-0, 


,13 


0.55 


26.21 


51609- 


■0304- 


■219 


SDSSJ141723.88-002220.8 


214, 


.349503 


-0.372432 


343, 


.356750 


55, 


,602268 


16.64 


1, 


.05 


-0, 


,31 


0.70 


23.05 


51613- 


■0305- 


■243 


SDSSJ142408.24-002930.7 


216, 


.034348 


-0.491848 


345, 


.608093 


54. 


,511181 


15.71 


1, 


,17 


-0, 


,02 


0.60 


25.01 


51613- 


■0305- 


■488 


SDSSJ142826.28-I-002915.4 


217, 


.109497 


0.487611 


348, 


.109955 


54. 


,617874 


17.54 


1, 


,17 


-0, 


,26 


0.60 


24.71 


51614- 


■0281- 


■438 


SDSSJ112513.99-I-004207.5 


171, 


.308289 


0.702082 


261, 


.294983 


56, 


,438610 


18.51 


1, 


.13 


-0, 


,23 


0.60 


24.60 



Note. — The first two columns are object names and the next four columns contains the astrometry (ra, dec, 1, b) for each object. The magnituc 
the Balmer lines. The positions are listed in the next five columns. The radial velocities and errors are listed next. The next three columns are a 
method and the scale width vs. shape method. The complete version of this table is in the electronic edition of the Journal. The printed edition ( 
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Table 3. Estimates of the Circular Velocities as a Function of Distance 



r (kpc) Vcir,l (km s ^) Vcir.ll (km s ^) 



7.5 


202 


± 


20 


215 


± 


20 


12.5 


227 


± 


20 


226 


± 


20 


17.5 


206 


± 


20 


180 


± 


20 


22.5 


170 


± 


20 


164 


± 


20 


27.5 


168 


± 


24 


183 


± 


20 


32.5 


162 


± 


27 


143 


± 


22 


37.5 


175 


± 


24 


183 


± 


39 


42.5 


207 


± 


30 


203 


± 


35 


47.5 


148 


± 


31 


166 


± 


30 


55.0 


180 


± 


54 


180 


± 


35 



Note. — Vcir,l denote the estimates based on 
Simulation I, while Y^ir II denote the estimates 
based on Simulation II. These estimates have been 
derived from the velocity scaling of the observed 
distribution that yields the best agreement to the 
simulated velocity distribution, as determined by a 
K-S test probability, and have been corrected for 
the effects of different halo stellar densities in the 
observations and simulations. Sec §3.2 for further 
information. 



